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ABSTRACT 

The Kepler Mission has acquired 33.5 d of continuous one- minute photometry of 
KPD 1946+4340, a short-period binary system that consists of a subdwarf B star 
(sdB) and a white dwarf. In the light curve, eclipses are clearly seen, with the deepest 
occurring when the compact white dwarf crosses the disc of the sdB (0.4 per cent) 
and the more shallow ones (0.1 per cent) when the sdB eclipses the white dwarf. As 
expected, the sdB is deformed by the gravitational field of the white dwarf, which 
produces an ellipsoidal modulation of the light curve. Spectacularly, a very strong 
Doppler beaming (also known as Doppler boosting) effect is also clearly evident at the 
0.1 per cent level. This originates from the sdB's orbital velocity, which we measure 
to be 164.0 ± 1.9 km s -1 from supporting spectroscopy. We present light curve models 
that account for all these effects, as well as gravitational lensing, which decreases the 
apparent radius of the white dwarf by about 6 per cent when it eclipses the sdB. We 
derive system parameters and uncertainties from the light curve using Markov Chain 
Monte Carlo simulations. Adopting a theoretical white dwarf mass-radius relation, the 
mass of the subdwarf is found to be 0.47 ± 0.03 M and the mass of the white dwarf 
0.59 ± 0.02 M . The effective temperature of the white dwarf is 15 900 ± 300 K. With 
a spectroscopic effective temperature of T e fj = 34 730 ± 250 K and a surface gravity of 
logg = 5.43 ± 0.04, the subdwarf has most likely exhausted its core helium, and is in 
a shell He burning stage. 

The detection of Doppler beaming in Kepler light curves potentially allows one to 
measure radial velocities without the need of spectroscopic data. For the first time, 
a photometrically observed Doppler beaming amplitude is compared to a spectro- 
scopically established value. The sdB's radial velocity amplitude derived from the 
photometry (168±4kms _1 ) is in perfect agreement with the spectroscopic value. Af- 
ter subtracting our best model for the orbital effects, we searched the residuals for 
stellar oscillations but did not find any significant pulsation frequencies. 

Key words: binaries: close - binaries: eclipsing - stars: individual (KPD 1946+4340) 
- subdwarfs - Kepler. 



1 INTRODUCTION 

Subdwarf B stars are mostly assumed to be extreme horizon- 
tal branch stars, i.e., cor e helium burning stars with a thin 
inert hydrogen envelope (|Hebe r 1986; S affer et alJ[l993 ). In 
order to reach such high temperatures and surface gravities, 
the progenitor must have lost almost its entire hydrogen en- 
velope. The majority of sdBs is expected to have lost its 
envelope v i a binary in teraction channels, as elaborated by 
lHan et al.1 (|200l l2003h . Our target, KPD 1946+4340 (KIC 
7975824), is a subdwarf B star (sdB) with a white dwarf 
(WD) companion in a 0.403739(8) day orbit (Morales- Rueda 
et al. 2003), which identifies the theoretical formation chan- 
nel for this system as the second common-envelope ejection 
channel of Han et al. (2002, 2003). In this scenario the white 
dwarf is engulfed by the sdB progenitor as it ascends the 
first giant branch. The white dwarf will deposit its angular 
momentum in the atmosphere of the giant and spin up the 
envelope until it is ejected. There are two subchannels to this 
scenario, depending on the initial mass of the progenitor. If 
sufficiently massive, it will ignite helium non-degeneratively, 
and the resulting extended horizontal branch (EHB) star 
will have a mass of ~ 0.35 M©. The more common scenario, 
starting with a roughly solar-mass giant, produces an EHB 
star with a mass that must be very close to the helium flash 
mass of 0.47 M©. A third possibility occurs when the white 
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dwarf companion ejects the envelope before the core has 
attained sufficient mass to ignite helium. In this case the re- 
maining core will evolve directly to the white dwarf cooling 
track. On its way it crosses the domain of the EHB stars, 
but without helium ignition the period in which it appears 
as an sdB star is brief, making this channel a very small 
contributor to the sdB population. For a recent extensive 
review on hot subdwarf sta rs, their evolution and observed 
properties, see iHeberl ([20091 ). 

The exact physical details involved in common-envelope 
ejection are not well understood. This uncertainty is com- 
monly embodied in the efficiency parameter a, which de- 
notes the amount of orbital energy used to eject the envelope 
(see e.g. Ide Koollll990l : iHu et alJl2Q07l ). Eclipsing subdwarf 
binaries could help constrain the permitted values of a, but 
studies have hitherto been hampered by the fact that both 
sdB+WD and sdB+M-dwarf binaries have virtually invis- 
ible companions and are therefore single lined, leaving the 
masses indeterminate. Firmly establishing the parameters of 
both components of a post-CE system therefore has substan- 
tial implications not just for confirming that our formation 
scenarios are correct, but also in order to tune future binary 
population synthesis studies by confining the a parameter. 

The target studied here, KPD 1946+4340, is an sdB 
star discovered by the Kitt Peak Downes survey ([Downed 
I1986T ). KPD 1946+4340 h as a V-band magnitude of 14.284 + 
0.027 (lAllard et al.lll994h . a ^/-magnitude of 14.299 ± 0.002 
( Wes emael et al.l |_1992) and a Kp (Kepler) magnitude of 
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14 .655. The star was in c luded in the radial velocity survey 
of iMorales-Rueda et al.l ([2003h , who found the target to be 
a spectroscopic binary with a period of 0.403739(8) d and 
a velocity amplitude K\ — 167 =b 2kms -1 . They also con- 
cluded that the sdB primary should be in a post-EHB stage 
of evolution, due to its relatively low surface gravity placing 
it above the canonical EHB in the HR diagram. This implies 
that the sdB exhausted all available helium in its core and 
is now in a shell helium burning stage. Assuming the sdB 
mass to be 0.5 M©, they found a minimum mass of 0.628 M© 
for the companion. 

In this paper we present the first light curve of 
KPD 1946+4340 obtained from space. The target was ob- 
served for 33.5 d by the Kepler Mission, and the light curve 
reveals sufficient low level features to permit purely photo- 
metric measurements of velocities, radii and masses of both 
components. A r eview of the Kepl er Mission and its first 
results is given in iKoch et al.1 (|2010h . 

We combine the Kepler photometry with new and old 
spectroscopic measurements, use light curve modelling to 
estimate the system parameters and Markov Chain Monte 
Carlo (MCMC) simulations to establish the uncertainties. 
The relativistic Doppler beaming effect is clearly detected 
in the light curve, and can be used to determine the orbital 
velocity of the primary. This effect, which is also known as 
Doppler boosting, was recently noted in a Kepler light curve 
of KOI- 74 by van Kerkwijk et al. (2010). We present the first 
comparison of a radial velocity amplitude as derived from 
the amplitude of the Doppler beaming to the spectroscopi- 
cally determined value. We also use the spectroscopic data 
to provide a revised ephemeris, as well as to determine the 
effective temperature, surface gravity and helium fraction 
of the atmosphere. After detrending the Kepler light curve 
with our best model for the orbital effects, we search the 
residuals for stellar oscillations. 



2 OBSERVATIONS, RADIAL VELOCITIES 
AND UPDATED EPHEMERIS 

We used 33.5 d of Ql short cadence Kepler data with a 
time resolution of 59 s. A review of the chara cteristics of the 
first sh ort cadence datasets is presented by iGilliland et al.l 
(l2010bh . The data were delivered to us through the KASOC 
(Kepler Asteroseismic Science Operations Center) websit^E 
The level of contamination of the fluxes by other stars is 
poorly known. We used the raw fluxes and assumed zero con- 
tamination, which is justified by the absence of other signif- 
icantly bright sources within ten arcsec of KPD 1946+4340. 
We applied a barycentric correction to the Kepler timings 
and converted them from UTC to barycentric dynamical 
time (TDB). The raw data show a ~ 2 per cent downward 
trend over the 33.5 d which we assume is instrumental in ori- 
gin. We removed this variation by fitting and dividing out 
a spline function. Out of the original dataset of 49170, we 
rejected 54 points because of a bad quality flag. Then, after 
initial light curve model fits to be described below, we re- 
jected another 87 points because they differed by more than 
3.5 cr from our model. The full light curve we used for our 

1 http://kasoc.phys.au.dk/kasoc/ 



analysis is shown on Fig. [TJ The time span of the dataset is 
BMJD(TDBj3 54 964.00314 to 54 997.49381. 

The high signal-to-noise spectra from I Green et al.l 

(2008) were used to derive T e ff and logg; see Section [H for 
details. 

To refine the orbi tal period determination of 
IMorales-Rueda et al ] (|2003l ). spectra were collected with the 
2.56-m Nordic Optical Telescope (NOT) on 2009 December 
5, 9 and 10. Using the ALFOSC spectrograph and a 0.5 
arcsec slit, we obtained eleven spectra with exposure times 
of 300 s and a of resolution R ~ 2000, covering a wavelength 
region of 3500 — 5060A. We measured the radial velo cities 
using mult i- Gaussian fits (|Morales-Rueda et~ al. 2003); the 
values are listed in Table [1] We first fitted the data of 
IMorales-Rueda et al.1 (|2003T ) and the new data separately to 
determine their RMS scatter. From these fits we found that 
it was necessary to add 3.6 and 4.1 kms -1 in quadrature to 
the uncertainties of the two datasets to deliver a unit \ 2 P er 
degree of freedom; these values probably reflect systematic 
errors due to incomplete filling of the slit. We then carried 
out a least-squares sinusoidal fit to the combined dataset 
finding a best fit of \ 2 = 21.6 (25 points). The next best 
alias had \ 2 — 39, and so we consider our best alias to 
be the correct one. This gave the following spectroscopic 
ephemeris: 

BMJD(TDB) = 53 652.84813(62) +0.40375026(16)^, 

marking the times when the sdB is closest to Earth. The 
corresponding radial velocity amplitude was K\ — 164.0 ± 
1.9 kms -1 . The WD's spectrum will slightly reduce the ob- 
served velocity amplitude of the sdB. We estimate that this 
effect is less than lkms -1 . The radial velocity measure- 
ments and the fit are shown in Fig. [2] From fitting light 
curve models to the Kepler photometry (see Section 13.4)) we 
derived the following photometric ephemeris: 

BMJD(TDB) = 54 979.975296(25) + 0.40375000(96)^. 

The two independent periods agree to within their uncer- 
tainties. The 9yr baseline of the spectroscopic ephemeris 
gives a more precise value and henceforth we fix it at this 
value in our light curve models. The cycle count between the 
spectroscopic and photometric ephemerides (using the spec- 
troscopic period) is 3287.0001 ± 0.0015, an integer to within 
the uncertainties. The Kepler-based zeropoint is the more 
precise one and is therefore retained as a free parameter in 
our models. 



3 LIGHT CURVE ANALYSIS 

The Kepler light curve we analyse in this paper (Fig. [TJ re- 
veals that KPD 1946+4340 is an eclipsing binary. We graph- 
ically determined the eclipse depths and durations. The 
eclipses of the WD by the sdB are 0.13 =b 0.03 per cent deep 
and the eclipses of the sdB by the white dwarf 0.38 ±0.03 per 
cent. The duration of the eclipses at half maximum depth is 
0.0236 ±0.0003 in orbital phase units. There is a clear asym- 
metric ellipsoidal modulation pattern in which the flux max- 
imum after the deeper eclipses is larger than the maximum 

2 BMJD(TDB) refers to Barycentric-corrected Modified Julian 
Date on the Barycentric Dynamical Timescale. 
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Table 1. Radial velocities of the sdB in KPD 1946+4340 deter- 
mined from NOT spectroscopy. 



BMJD (TDB) 


RV (kms 


- 1 ) 


55170.81118 


-139.8 


± 


3.5 


55170.86268 


-162.0 


± 


3.3 


55170.90743 


-109.2 


± 


3.8 


55170.95896 


22.5 


± 


5.7 


55174.80952 


-73.5 


± 


3.0 


55174.84999 


-154.7 


± 


3.3 


55174.90291 


-171.9 


± 


4.0 


55174.93220 


-127.6 


± 


4.7 


55175.80006 


2.3 


± 


4.4 


55175.84379 


102.4 


± 


4.7 


55175.88328 


148.3 


± 


6.5 



after the shallower eclipses. We attribute this to Doppler 
beaming, see Section [3721 

To determine the system properties, we modelled the 
light curve with the LCURVE code written by TRM ( for a 
description of the code, see Co pperwheat et al. I Hold, Ap- 
pendix A). This code uses grids of points to model the two 
stars, taking into account limb darkening, gravity darkening, 
Doppler beaming and gravitational lensing when the white 
dwarf eclipses its companion. It assumes the ellipsoidally de- 
formed star to be in corotation with the binary orbit, which 
is usually a good assumption because of the large tidal inter- 
actions between the two binary components. Reprocessing 
of light from the sdB by the WD is included in the light 
curve models as well ("reflection effect"). 

To speed up the computation of the models used in this 
paper we implemented a new option whereby a finer grid 
can be used along the track of the white dwarf as it eclipses 
the sdB. This reduces the overall number of points needed to 
model the light curve to the demanding precision required 
to model the Kepler data. In addition, we only used the 
finely-spaced grid during the eclipse phases, taking care to 
make the model values continuous when changing between 



200 




0.0 0.2 0.4 0.6 0.8 1.0 
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Figure 2. Radial velocity curve of the sdB in KPD 1946+4340. 
Both our new radial veloci ty measurements and the ones from 
iMorales-Rueda et al.l (120031 ) are shown, folded on the orbital pe- 
riod. The error bars of the datapoints show the uncertainties after 
adding 3.6 kms -1 an d 4.1 k ms -1 in quadrature to the values of 
IMora les-Rued a et al.l (|2003l ) and our new data respectively. We 
find a radial velocity amplitude of K\ = 164.0 ± 1.9 kms -1 . 

grids by applying normalisation factors (very close to unity) 
to the coarse grid fluxes. We used ~ 100 000 (~ 37000) grid 
points for the fine (coarse) grids for the sdB and 3000 for the 
white dwarf. To model the finite exposures more accurately 
during the eclipses, where smearing occurs due to the 1 m 
integration time, we calculated 7 points for each exposure 
(i.e., one point for every ~ 10s) and took their trapezoidal 
average. 

3.1 Gravity darkening and limb darkening 
coefficients 

We used model spectra to compute the gravity darkening 
coefficient (GDC) of the sdB and the limb darkening coeffi- 
cients for both the sdB and the white dwarf, which are all im- 
portant parameters for the modeling of a close binary's light 
curve. The GDC is needed to model the effects of the white 
dwarf's gravity on the sdB, which slightly distorts the sdB's 
shape. The bolometric flux from a stellar surface depends 
on the local gravity as T 4 oc g^ b in which /3b is the bolo- 
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metric GDC. For radiative stars, /3b = 1 (jvon Zeipellll924h . 
We observe the band- limited stellar flux, not the bolomet- 
ric flux, and hence we require a different coefficient defined 
by / oc g^ K . The GDC for the Kepler bandpass, /3k was 
computed from 

d log / 



Fx 



1 - B 



t) 



(2) 



_ dlogl | a log J dlogT 
dlogg d log g dlogT dlogg 



(1) 



in which / is the photon-weighted bandpass-integrated spe- 
cific intensity at fi = 1 and ^2|Z = = 0.25. We used 
a grid of sdB atmosphere models calculated from the LTE 
model atmosp here grid of [ Heber et al.l (|2000h using the Lin- 
for program (|Lemkd \l9m and assumed T e ff = 34 500K, 
\ogg = 5.5, log(n He /n H ) = -1.5 and log (Z/Z®) = -2. 
To estimate the interstellar reddening, w e compared the ob - 
served B — V colour of -0.20 ±0.01 mag (|Allard et al.lll994T ) 
with the colours expected from a model atmosphere. We 
found an intrinsic colour of B — V = —0.26 mag and conse- 
quently adopted a reddening of E(B — V) = 0.06. To account 
for this interst ellar reddening the m odel spectra were red- 
dened following ICardelli et afl(jl989h . The gravity darkening 
coefficient was found to be /3k = 0.448. 

Using a model for the same set of parameters, we com- 
puted limb darkening coefficients for the sd B. We a d opted 
the 4-parameter limb darkening relation of Claret (2004, 
equation 5) and determined a± = 0.818, a2 = —0.908, 
a 3 = 0.755 and a 4 = -0.252. 

For the white dwarf, angle- dependent model spe ctra 
were calculated using the code of iGansicke et al.l (|l995h for 
T e ff = 17 000 K (estimated from a comparison of model sur- 
face brightnesses given initial light curve fits) and logg = 
7.8. We adopted the same limb darkening law as for the 
sdB and found ai = 0.832, a 2 = -0.681, a 3 = 0.621 and 
a 4 = -0.239. 



3.2 Doppler beaming factor 

The asymmetry in KPD 1946+4340's ellipsoidal modulation 
pattern is the result of Doppler beaming. Doppler beaming 
is caused by the stars' radial velocity shifting the spectrum, 
modulating the photon emission rate and beaming the pho- 
tons somewhat in the direction of motio n. The effect was, as 
far as we are aware, first discussed in iHills fe Dalel dl974h 
for ro tation of white dwarfs and by IShakura & Postnovl 
(Il987h_ for orbital motion in binaries. It was first observed 
by iMaxted et all fcOOOh . Its expected dete ction in Kepler 
light c urve s was suggested and discussed bv lLoeb & Gaudi 
<|20Q3h and lZucker et al.1 <|2007h . Van Kerkwiik et al. fcOlOh 
report the detection of Doppler beaming in the long ca- 
dence Kepler light curve of the binary KOI- 74. For the first 
time, they measured the radial velocity of a binary compo- 
nent from the photometrically detected beaming effect. The 
measured radial velocity amplitude, however, did not match 
the amplitude as expected from the mass ratio derived from 
the ellipsoidal modulation in the light curve. The derived 
velocity of the primary of KOI- 74 is yet to be confirmed 
spectroscopically. For KPD 1946+4340, radial velocities are 
available which allows the first spectroscopic check of a pho- 
tometrically determined radial velocity. 

For radial velocities that are much smaller than the 
speed of light, the observed flux F\ is related to the emitted 
flux Fq \ as 



with B the beami ng factor B — 5 + din Fa /din A 
([Loeb &; Gaudill2003h . The beaming factor thus depends on 
the spectrum of the star and the wavelength of the obser- 
vations. For the broadband Kepler photometry, we use a 
photon weighted bandpass-integrated beaming factor 



(B) 



j exXFxBdX 
J exXF x dX 



(3) 



in which e\ is the response function of the Kepler bandpass. 

We determined the beaming factor from a series of fully 
metal line-blanketed LTE models ([Heber et al hood , see also 
Section [3. 1|) with metallicities ranging from log (Z/Z@) = 
—2 to +1, as well as from NLTE models with zero metals and 
with Blanchette metal composition (see Section [4] of this pa- 
per for more information about the NLTE models) . Without 
taking reddening into account, the beaming factor is found 
to be (B) = 1.30 ± 0.03. The uncertainty incorporates the 
dependence of the beaming factor on the model grids and 
the uncertainty on the sdB's effective temperature, gravity 
and, most importantly, metallicity. The metal composition 
of the model atmospheres is a poorly known factor that can 
only be constrained with high-resolution spectroscopy. 

This time, the effect of reddening has to be accounted 
for by changing the spectral response accordingly instead 
of reddening the model atmosphere spectrum. Using a red- 
dened spectrum would in this case erroneously imply that 
the reddening is caused by material that is co-moving with 
the sdB star. With reddening, the beaming factor is deter- 
mined to be ~ 0.006 lower. Reddening thus only marginally 
affects the beaming of KPD 1946+4340 but should certainly 
be taken into account in case of higher reddening values. 

There are three contributions to the beaming factor. 
The enhanced photon arrival rate of an approaching source 
contributes +1 to the beaming factor. Aberration also in- 
creases the number of photons that is observed from an 
approaching source, adding +2 to the beaming factor be- 
cause of the squared relation between normal angle and 
solid angle. Finally, when the sdB comes towards us, an 
observed wavelength A corresponds to an emitted wave- 
length A e = A (1 -\-v r /c). Since sdBs are blue, looking at 
a longer wavelength reduces the observed flux which coun- 
teracts the other beaming factor components. In case of an 
infinite temperature Ray leigh- Jeans spectrum this Doppler 
shift contribution to the beaming factor would be —2. For 
the primary of KPD 1946+4340, we find a contribution of 
~ —1.70 which brings the total beaming factor to ~ 1.30. 
The contribution of the Doppler shift does not always have 
to be negative; a red spectrum could actually increase the 
effect of beaming. 



3.3 Light curve model 

A typical fit to the data is shown in Fig. [3] with the different 
contributions switched on, one- by-one. From the residuals 
(bottom panel) it is clear that the model reproduces the 
variations at the orbital period very well. When we fit the 
light curve outside the eclipses with sine curves to represent 
the reflection effect, ellipsoidal modulation and the beaming, 
the phase of the ellipsoidal modulation is found to be off 
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by 0.0072 ±0.0010 in orbital phase units. We do not know 
the origin of this offset, which also gives rise to the shallow 
structure that is left in the residuals. Our best fits have 
X 2 = 52032 for 48929 data points. 

The significance of the Doppler beaming is obvious, and 
even the more subtle gravitational lensing effect is very sig- 
nificant, although it cannot be independently deduced from 
the data since it is highly degenerate with changes in the 
white dwarf's radius and temperature. One part of the grav- 
itational lensing is caused by light from the sdB that is 
bent around the white dwarf, effectively making the white 
dwarf appear smaller. The second lensing contribution is a 
magnification effect which is caused by the altered area of 
the sdB that is visible, given that surface brightness is con- 
served by the lensing effect. In the case of KPD 1946+4340, 
the first part of the lensing is the most important. Lensing 
effects in compact b i naries were d iscus sed in e.g. Maederl 
(Il973h, iGouldl (Il995h. IMarshl (|200lh and|Agoj (|200ll2003l ). 
Sahu & GillilanJ (2003|) explored the expected influence of 
microlensing effects on light curves of compact binaries and 
planetary systems observed by Kepler. They found that the 
lensing effect of a typical white dwarf at 1 AU of a main 
sequence star will swamp the eclipse signal. A transit of a 
planet, which is of similar size but a lot less massive, can 
therefore easily be distinguished from an eclipse by a white 
dwarf. In the case of KPD 1946+4340, the separation of the 
two components is a lot less. An eclipse is still seen, but with 
reduced depth. For the most likely system parameters, grav- 
itational lensing reduces the eclipse depth by ~ 12 per cent 
which is equivalent to a ~ 6 per cent reduction of the ap- 
parent white dwarf radius. The effect of gravitational lensing 
is implement ed in our light curve modelling code following 
IMarshl (l200lh . 



3.4 Markov Chain Monte Carlo simulation 

The parameters which determine models can be fixed by 
minimisation of % 2 - If the signal-to-noise is high, a quadratic 
approximation around the point of minimum % 2 can lead to 
the uncertainties of, and correlations between, the best-fit 
parameters. The Kepler data have superb signal-to-noise, 
but owing to the very shallow depths of the eclipses the 
quadratic approximation does not work well. Strong corre- 
lations between several parameters play a significant role 
in this problem. The duration of the eclipses essentially fix 
the scaled radius of the sdB star (which we take to be the 
primary) n = Ri/a, where a is the binary separation. The 
scaled radius is a function of orbital inclination i, r± = r\(i). 
The depth of the eclipse of the sdB by the white dwarf fixes 
the ratio of radii rijr\ = R2/R1, so r2 is also a function of 
orbital inclination. The duration of the ingress and egress 
features provides an independent constraint on r2 as a func- 
tion of i, which can break the degeneracy. In this case, how- 
ever, one is limited by a combination of signal-to-noise and 
the minute-long cadence which is not sufficient to resolve 
the ingress/egress features. 

Under these circumstances, a Markov Chain Monte 
Carlo (MCMC) method can be very valuable. The MCMC 
method allows one to build up a sequence of models in which 
the fitting parameters, which we denote by the vector a, 
have a probability distribution matching the Bayesian pos- 
terior probability of the parameters given the data, P(a|d). 




1.0 

Orbital phase 

Figure 3. Phase-folded light curve of KPD 1946+4340 (green, 
datapoints grouped by 30) and our best fitting model (black). In 
the top panel, only the eclipses and reflection effects are modelled. 
In the second panel, ellipsoidal modulation is added. In the third 
panel, gravitational lensing is taken into account as well, which 
affects the depth of the eclipse at orbital phase 0.5. The bottom 
panels show the full model - taking into account Doppler beaming 
- and the residuals (grouped by 30 in green and grouped by 600 
in black). 



From long chains of models one can then calculate variances 
and plot confidence regions. The MCMC method also has 
the side benefit of helping with the minimisation which can 
become difficult when parameters are highly correlated. For 
data in the form of independent Gaussian random variables, 
this probability can be written as 



P(a|d) oc P(e 



-xV2 



(4) 



i.e., the product of one's prior knowledge of the model pa- 
rameters and a factor depending upon the goodness of fit as 
expressed in % 2 . We implemented t he MCMC method fol- 
lowin g procedures along the line of Collier Cameron et al. 
(J2007|). We incorporated prior information in two ways. In 
all cases we used our constraint K\ — 164 ± 2kms _1 . 
Using our own spec t roscop ic analysis and the results of 
iMorales-Rueda et al.l ([2003) , we decided to put also the fol- 
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lowing constraint on the effective temperature of the sdB: 
Ti = 34 500 ± 400 K. These two constraints were applied by 
computing the following modified version of \ 2 

oi (t>( \a\\ 2 , ( K x -\^ \ 2 ( T\ — 34 500 \ 2 
-21n(P(a|d)) = X — J + ( 400 J ' (5) 

where K\ and Ti are the values in the current MCMC model 
under test. The period of the binary orbit was kept fixed 
at the spectroscopically determined value (see Section [2]). 
The parameters that were kept free during the modelling 
are the scaled stellar radii Ri/a and it^/a, the mass ratio 
q, the inclination i, the effective temperature of the WD 
T2, the beaming factor and the zero point of the ephemeris. 
The radial velocity scale (which leads to the masses M\ and 
M2) and the effective temperature of the primary Ti were 
included in the fits as well, but with the spectroscopically 
allowed range as a prior constraint. Note that the beaming 
factor (B) is kept as a free parameter, which allows the code 
to fit the Doppler beaming amplitude while we constrain the 
allowed range of K\. By comparing the MCMC results for 
(B) with the theoretical beaming factor, we can check if the 
beaming amplitude is consistent with our expectations. 

As explained above, this led to parameter distributions 
with strong correlations between R±, R2, Mi, M2, q, etc. The 
mass-radius relations for the two stars are shown in Fig. [4l 
The favoured parameters for the secondary (left panel) 
clearly show that it is a white dwarf and the allowed distri- 
bution nicely crosses the expected track of mass-radius (solid 
green line) which we calculat ed from the zero temperatur e 
relation of Eggleton (quoted in lVerbunt &; Rappaportf lQSS), 
inflated by a factor o f 1.08. We estimated this f actor from 
the cooling models o f iHolberg fc Bergeron! (I2006T ) for white 
dwarfs with a mass between 0.5 and 0.7 M (white dwarfs 
with masses outside this range are ruled out by the mass- 
radius relation) using our preferred temperature for the 
white dwarf of around 16 000 K and assuming an envelope 
that consists of M H = 10~ 4 M W d and M He = 10~ 2 M W d- 

Given that the secondary is a white dwarf and given 
that the secondary's theoretical mass-radius relation inter- 
sects the mass-radius distribution very nearly within the 1-a 
region, we also undertook MCMC runs where we added the 
prior constraint that the secondary had to match the white 
dwarf M-R relation to within an RMS of 5 per cent. This 
was added in exactly the same manner as the K\ and Ti 
constraints. The system parameters we derive from these 
MCMC runs are listed in Table [2] The mass-radius rela- 
tion after applying the constraint is shown in Fig. [5] Espe- 
cially after applying the white dwarf mass-radius relation 
constraint, the sdB's mass-radius relation fits perfectly with 
the one defined by the surface gravity derived from spec- 
troscopy in Section [H The correlation coefficients between 
the different parameters are given in Table [3j The binary's 
inclination, its mass ratio and the stellar radii are highly 
correlated. 

3.5 Variability in residuals 

One of the goals of the Kepler Mission is to allow detailed 
asteroseismic studies of pulsati ng stars. The a s terosei smol- 
ogy programme is discussed in iGilliland etHI (|2010aT ). For 
more information about the search for pulsation s in compact 
objects with Kepler, see l0stensen et al.l (|2010al ). 



Table 2. Properties of KPD 1946+4340. The orbital period and 
the effective temperature of the sdB were derived from spec- 
troscopy. The other parameters are obtained by modelling the 
Kepler light curve. The uncertainties on these values are deter- 
mined by MCMC analysis, using the prior constraint that the 
white dwarf mass-radius relation has to match the Eggleton rela- 
tion to within 5 per cent RMS. 



Primary (sdB) Secondary (WD) 



^orb (d) 


0.40375026(16) 


q 


1.27 + 0.06 


i (deg) 


87.14 + 0.15 


R(Rq) 


0.212 ± 0.006 0.0137 ± 0.0004 


M (M ) 


0.47 + 0.03 0.59 + 0.02 


T eff (K) 


34 500 + 400 15 900 + 300 



Table 3. Correlation coefficients of the different parameters that 
were varied in the MCMC simulations, after applying the Eggle- 
ton mass-radius relation constraint. 





R2 


i 


Ti 


T 2 


Q 


Ri 


0.95 


-0.95 


0.02 


0.02 


-0.95 


R2 




-0.98 


0.07 


0.02 


-0.99 


i 






-0.06 


-0.02 


0.96 


Ti 








0.44 


-0.09 


T 2 










-0.02 



The Fourier transform of the original light curve of an 
eclipsing close binary like KPD 1946+4340 is highly contam- 
inated by frequencies and their harmonics due to the binary 
orbit. Subtraction of a good model of the binary signatures 
of the light curve allows one to get rid of this contamina- 
tion. Since a number of sdBs have been found to be multi- 
periodically pu lsating (for a r eview on asteroseismology of 
EHB stars, see l0stensenll2OO9h . we checked the residuals of 
the light curve for signs of pulsations. Usi ng the analysis 
metho d and significance criteria outlined in lDegroote et aD 
(2009), 8 significant frequencies were found, which are listed 
in Table H 

fi is a known artefact frequency caused by an eclips- 

Table 4. Significant variability frequencies in the residuals of the 
light curve of KPD 1946+4340. The signal to noise (S/N) value 
was determined by dividing the amplitude of the peak by the 
uncertainty on the amplitude. 





Frequency (d 1 ) 


Amplitude 
(/xmag) 


S/N 




A 


0.2758 + 0.0011 


97.5 + 7.2 


13.5 


= /2/2 


h 


0.5936 + 0.0013 


97.5 + 7.2 


13.5 


instrumental 


h 


0.1417 + 0.0014 


86.6 + 7.2 


12.1 


= /a/4 


h 


1.1820 + 0.0015 


81.0 + 7.2 


11.3 


= 2/2 


h 


1.7730 + 0.0018 


67.5 + 7.2 


9.4 


= 3/2 


h 


440.4386 ± 0.0022 


54.2 + 7.2 


7.6 


instrumental 


h 


4.9548 ± 0.0024 


49.1 + 7.2 


6.9 


2 J ^ orb 


h 


0.3115 + 0.0027 


44.9 + 7.2 


6.3 
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0.4 0.6 0.8 1.0 1.2 1.4 0.5 1.0 1.5 2.0 



M W D (M ) M sdB (M ) 

Figure 4. Mass-radius relations of the WD (left) and sdB (right). One tenth of our MCMC models are shown (black dots). The contour 
plots show the regions in which 68, 95 and 99 per cent of the models reside. The contours that are shown are somewhat artificially 
broadened by the binning process. The Eggleton mass-radius relation, inflated by a factor 1.08 to allow for the finite WD temperature 
(see text for details), is shown as a solid green track on the left plot. The mass-radius relation intersects the Eggleton relation very 
nearly within its 1-a region. On the right panel the solid green line gives the mass-radius relation for logg = 5.45, the dashed lines for 
\ogg = 5.40 (left) and \ogg = 5.50 (right). 




ing binary that was used as one of the fine- guidanc e stars 
during Ql fsee lHaas etaHl20ld : [jenkins et al]l2010h . Four 
other frequencies (/i, /3, f± and fe) are related to fc. The 
highest frequency, fe, is related to the pr ocessing of the long 
cadence data fsee iGilliland et aDfeOlObh . f*i is the first har- 
monic of the orbital frequency of KPD 1946+4340, which 
indicates that there is still a weak orbital component left 
after subtracting our light curve model, fs is not related to 
any of the other frequencies and corresponds to a period that 
is too long to arise from stellar pulsations of the WD or the 
sdB. If it is real, the signal might result from the rotation 
of the WD, or it might be due to a background star. 

The best candidate peak for p— mode pulsations of the 
sdB is at 5018.2 /iHz with an amplitude of 37 /xmag, but fur- 
ther data is needed to confirm that the sdB is pulsating. 



From ground based data, l0stensen et al.1 (| 201 Obi ) did not 
detect pulsations, with a limit of 0.68 mma. This is consis- 
tent with the Kepler photometry. 



4 SPECTROSCOPIC ANALYSIS 

Low resolution high S/N spectra for KPD 1946+4340 were 
taken with the B&C spectrograph at Steward Observatory's 
2.3-m Bok telescope on Kitt Peak, as par t of a long term ho - 
mogeneous survey of hot subdwarf stars ([Green et al.ll2008h , 
in September and October 2004. The spectrograph param- 
eters, observational procedures, and reduction techniques 
were kept the same for all the observing nights. 

The 400/mm grating, blazed at 4889 A, gives a resolu- 
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tion of R = 560 over the wavelength region 3620 — 6895 A, 
when used with the 2.5 arcsec slit. The spectra were taken 
during clear or mostly clear conditions with integration 
times between 1050 and 1200 s. Approximately 1000 bias and 
flat-field images were obtained for each run. The data reduc- 
tions were performed using standard IRAF tasks, and each 
night was flux-calibrated separately. The individual spectra 
were cross-correlated against a super-template to determine 
the relative velocity shifts, and then shifted and combined 
into a single spectrum. Although the resolution is rather low, 
the S/N is quite high: 221/pixel or 795/resolution element 
for the combined spectrum. The continuum fit to the com- 
bined flux- calibrated spectrum were done with great care to 
select regions devoid of any weak lines, including expected 
unresolved lines of heavier elements. 

The final KPD 1946+4340 spectrum was fitted using 
two separate grids of NLTE models designed for sdB stars, 
in order to derive the effective temperature, surface grav- 
ity and He/H ratio. The first set of models assumed zero 
metals, while the second included an adopted distribution 
of meta ls based on the analysis of FUSE spectra of five sdB 
stars b y felanchette et al.l (|2008h , see also lVan Grootel et al.l 
(|2010h . From the set of models without metals, we derive 
log# = 5.45 ± 0.04, T eff = 34 400 ± 220 K and log(He/H) = 

— 1.37 ±0.05. Assuming the Blanchette composition, we find 
log g = 5.43 ± 0.04, T eff = 34 730 ± 250 K and log(He/H) = 

— 1.36 =b 0.04. These results are in good agreement with 
log<? = 5.37 ± 0.10, T eff = 34 500 ± 1000 K, log(He/H ) = 
-1.35±0.10 determined bv lMorales-Rueda et alT(|2003[ ) and 
logfl 5 - 4 , 3 ± Q- lQi T eff = 34 200 ± 500 K determined by 
iGeier et al.1 (linoh using different spectra and model grids. 

The fit definitely improves when going from the zero- 
metal solution (Fig. [6} to the Blanchette composition 
(Fig. [7]), although there still remains a slight "Balmer" 
problem, especially noticeable in the core of H/3. There are 
definitely metals in the spectrum of KPD 1946+4340: the 
strongest features are 1) an unresolved C III + Nil complex 
around 4649 A (compare the two figures for that feature), 
and 2) another weaker complex (C III + O II) in the blue 
wing of H5 that the Blanchette model reproduces quite well. 
All of the major discrepancies between the spectra and the 
models are due to strong interstellar absorption: the K line 
of Ca II in the blue wing of He, the Ca II H line in the core of 
He, and the Nal doublet strongly affecting the red wing of 
He I 5876. It is reassuring that the derived atmospheric pa- 
rameters are not too strongly dependent on the presence of 
metals, as might be expected for such a hot star, particularly 
one in which downwards diffusion of metals is important. 



5 DISCUSSION 

The beaming factor we derived for KPD 1946+4340 using 
MCMC runs is (B) = 1.33 ± 0.02, which is in perfect agree- 
ment with the theoretically expected value calculated in Sec- 
tion [32] The uncertainty on the beaming factor is a direct 
reflection of the uncertainty on the spectroscopic radial ve- 
locity amplitude of the sdB. If, contrary to our assumption, 
the Kepler fluxes would be severely contaminated by light 
from other (constant) stars, the observed beaming factor 
would be lower. The distribution of beaming factors from 
our MCMC computations is shown in Fig. [8] If the radial ve- 



25(H) 




1.25 1.30 1.35 1.40 

Beaming factor 



Figure 8. Distribution of the sdB's beaming factor for an MCMC 
run (using the M-R constraint for the WD; see text for details). 
The beaming factor is found to be (B) = 1.33 ± 0.02, which is in 
agreement with the theoretically expected (B) = 1.30 ± 0.03. 

locity would be measured from the Doppler beaming ampli- 
tude, using the theoretical beaming factor, the result would 
be 168 ± 4kms -1 compared to 164.0 ± 1.9 km s -1 derived 
from spectroscopy. The uncertainty on the photometric ra- 
dial velocity is dominated by the uncertainty on the theoret- 
ical beaming factor, primarily due to its dependence on the 
poorly known metallicity of the sdB, and to a lesser extent 
due to the uncertainties on the sdB's effective temperature 
and surface gravity. 

Under the assumption of corotation, we find a projected 
rotational velocity of the sdB of i>sin(z) = 26.6 + 0.8 kins" 1 . 
From spectrosco py and using LTE models with ten times 
Solar metallicity, IGeier et~al ] (|2010h found v sin(i) = 26.0 ± 
1.0 km s -1 , which is in agreement with our photometric re- 
sult. We conclude that the assumption of corotation is likely 
to be correct. 

The spectroscopically determined surface gravity of the 
sdB (logp = 5.43 ± 0.04 and 5.45 ± 0.04 using atmosphere 
models with and without metals respectively) agrees per- 
fectly with the surface gravity of 5.452 ± 0.006 we derived 
from the mass-radius distribu tion of our light curve models . 

As concluded earlier by iMorales-Rueda et al.1 (|2003h , 
the sdB is probably in a post-EHB phase. This is illustrated 
in Fig. [9j which shows the zero age extended horizontal 
branch (ZAEHB) and the terminal age extended horizon- 
tal branch (TAEHB) for an sdB with a typical core mass 
of 0.47 M©, together with evolutionary tracks for different 
hydrogen envelope masses (10~ 4 , 10~ 3 , 2 x 10~ 3 , 3 x 10" 3 
and 4 x 1O" 3 M ) from lKawaler fe Hostler! (|2005l ). 

Because of its low surface gravity, the sdB component of 
KPD 1946 falls in a region of the T e ff — logg plane relatively 
far from the center of the instability strip. However, at least 
one pulsator exists in this region of the T e ff — lo gg plane, 
V338 Ser, that should be in a post-EHB phase (see 0st enseni 
120091 Fig. 3). Moreover, 'transient pulsators' with varying 
pulsation amplitudes that can go down to undetectable val- 
ues in a pa rticular epoch might e xist (see the case of KIC 
2991276 in l0stensen et alJ l2Q10ah . For these reasons, and 
because we found at least one candidate p— mode pulsation 
frequency, it is worth continuing a photometric monitoring 
by Kepler. 

The white dwarf mass implies that it is a CO white 
dwarf. The progenitor of the sdB must have been the less 
massive star in the original binary and by the time it reached 
the ZAEHB, the white dwarf was already cooling. The ac- 
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Figure 6. Fit (bold lines) to the spectral lines of KPD 1946+4340 using NLTE sdB models assuming zero metals. 
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Figure 7. Fit (bold lines) to the spectral lines of KPD 1946+4340 using NLTE sdB models assuming a Blanchette metal composition. 



cretion of material by the white dwarf does not change 
the WD's internal energy content sig nificantly (see e.g. re- 
lated work on cataclysmic variables bv iTownslev &; Bildstenl 
120021 ). The cooling time of the white dwarf therefore sets 
an upper limit to the time since the sdB was on the ZA- 
EHB. For our best estimates of the tem perature and mass of 
the w hite dwarf, the cooling tracks of iHolberg fc Bergeron! 
(2006) indicate that it has been cooling for about 155 to 
170 Myr (depending on the unknown envelope composition). 
The sdB's evolution from the Zero Age Extended Horizontal 
Branch to its current pos t-EHB-phase took 125 to 145 Myr 
([Kawaler &; Hostlerl 120051 . depending on the exact current 
evolutionary stage), which means that the sdB must have 
formed very shortly after the white dwarf. 



6 SUMMARY 

We have analysed a 33.5-d short cadence Kepler light curve 
of KPD 1946+4340, as well as low resolution spectroscopy. 
In the light curve, primary and seconday eclipses, ellipsoidal 
modulation and Doppler beaming are detected. We model 
the binary light curve, taking into account the Doppler 
beaming and gravitational lensing effects. System param- 
eters and uncertainties are determined using Markov Chain 
Monte Carlo simulations. 

The binary is found to consist of a 0.59 ±0.02 M white 
dwarf and a 0.47 ±0.03 M post-EHB sdB star. The surface 
gravity and corotation rotational velocity of the sdB as de- 
rived from the light curve models are found to be consistent 
with spectroscopic values. The observed Doppler beaming 



© 2010 RAS, MNRAS 000.(110 



Kepler observations of the beaming binary KPD 1946+4340 11 



4.6 



4.8 - 




50000 45000 40000 35000 30000 25000 20000 



T eS (K) 

Figure 9. The sdB of KPD 1946+4340 in the T eff - \ogg plane. 
The theoretical zero age and terminal age extended horizontal 
branches for a 0.47 Mq are shown, together with evolutionary 
tracks for different envelope thicknesses (10 -4 , 10 -3 , 2 x 10 -3 , 
3 x 10 -3 and 4 x 10 -3 M ). The sdB is found to be in the post- 
EHB phase. 

amplitude is in perfect agreement with the amplitude ex- 
pected from spectroscopic radial velocity measurements. It 
would thus have been possible to derive the radial velocity 
amplitude of the sdB from the Kepler light curve directly. 

Subtracting a good light curve model allowed us to 
search for stellar oscillations. No significant stellar variabil- 
ity of the sdB or white dwarf could be detected yet. At least 
one candidate p— mode pulsation frequency was found, how- 
ever, and the sdB can also possibly be a transient pulsator. 
KPD 1946+4340 continues to be observed by Kepler. 
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